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We investigate, through Monte-Carlo simulations, the nature of the second order point in a Z 2 
(Bosonic) -I- Z 2 gauge theory in four dimensions. Detailed analysis of the critical exponents point 
to the Ising universality class. Relevancy to extended models and possible Non-Gaussian behavior 
is discussed. 


In the present paper we investigate the character of the 
critical point C (Fig 1) and show that the segment AC is 
mapped on the segment H = 0, T G [0,Tc] of the usual 
Ising model. We find that the Ising critical exponents 
and the Ising universal function are reproduced to a very 
high precision. 


I. INTRODUCTION 


The Z 2 gauge system has been investigated using 
mean-field method Q , Monte Carlo simulations |^,|| and 
exact solutions HJstl. The objective of these studies was 


to find the phase diagram from the action (equation 1.1), 
and to classify the order of the phase transitions. 


A. The Model 

The action of the Z 2 gauge lattice can be written as 

S = —fJp ^2 UijUjkUkiUii - (3 l X (1-1) 

° (hi) 

where Uij is the value of the link connecting the sites i 
and j and cr^ is the value of site i. The first sum is over 
all sets of four neighboring links (plaquettes) and the 
second sum is over all nearest neighbor sites. Both Uij 
and (Ti can obtain the values ±1 with Uij = Uji. The 
free energy is then 

F=^lnZ ( 1 . 2 ) 

with the partition function Z given by 

Z= e-^ (1.3) 

Uij ,ai 

and N is the number of site s in the lattice. 

From the action Eq. ( |l.l[ ) we look at two gauge invari¬ 
ant observables L and P define as 

P = {UijUjkUkiUH)= (1.4) 

L= (UiUijaj) = (1.5) 


where the averaging is over all sites and links of the lat¬ 
tice. The factors | and 4 are the ratios between the 
number of sites to the number of links and plaquettes in 
a four-dimensional lattice. 

In this system, L and P can be used as order param¬ 
eters replacing the standard magnetization which is not 
gauge invariant. They are order parameters in the sense 
that they exhibit singularities of the bulk thermodynam¬ 
ics. However they lack the property of a magnetization 
in that they never vanish (except at temperature zero) 
or zero. 



FIG. 1. Schematic plot of the phase diagram. The dashed 
lines represent a first order phase transition, and the contin- 
nons line represents a second order transition. The critical 
point is marked with G and the triple point by T. 


The phase diagram is shown schematically in Fig. [^. 
The three transition lines experienced by L and P meet at 
the triple point(A) and split the phase space into three 
regions. This phase space includes two familiar limits. 
In the limit /3p ^ 00 all plaquettes are forced to one. 
The gauge degrees of freedom therefore disappear and 
the system reduces to that of a normal four-dimension 


1 








Ising model. Here L experiences a second order transi¬ 
tion while P remains smooth. It has been shown ^ that 
the second order line is of Ising type all the way to the 
triple point. 

On the /3l = 0 axis the action reduces to a pure gauge 
theory and P experiences a first order transition. This 
behavior continues all the way up to the triple point A 
with L remaining smooth [Q. A line of first order in both 
L and P is then formed which ends with the critical point 
C. 


II. NUMERICAL SIMULATION ON THE 
LATTICE 

We implemented a standard Monte Carlo procedure 
using a heat bath algorithm to generate the configura¬ 
tions of the lattice. In this paper we shall use the term 
iterationto represent an update of all links and sites in 
the lattice. 

Even using this method, the simulations require sig¬ 
nificant computational resources restricting us to modest 
lattice sizes for our investigation. Therefore we used lat¬ 
tice sizes ranging from 4^ sites ( which other studies have 
shown to give good results up to 14“^, and at the 

critical point up to 18^. In order to minimize surface 
effects a toroidal geometry was adopted. 

Each simulation consisted of measurements of the ob¬ 
servables along a line of points in the phase diagram. 
Many such lines were made in order to cover the param¬ 
eter space and measure variations along different direc¬ 
tions. At each point the system was allowed to evolve 
for at least 10,000 iterations, with mean values of the 
observables calculated after discarding the first 1000. In 
the vicinity of the critical point, the relaxation time and 
fluctuations increased dramatically requiring more itera¬ 
tions to give reliable results. In this region we generated 
up to 100,000 iterations. 

Along the first order line AC, the behavior of L is like 
a magnetization with two metastable states L'^ and L~. 
However, unlike the Ising model, these two states are not 
symmetrical about the origin. By evaluating the mid 
point Lq between L+ and L~. we can consider L — Lq 
which is then symmetrical about zero. However Lq can 
only be found along the line of first order where there are 
two stable states present. 

Direct measurements of L are further complicated by 
tunneling between L+ and L~. This has the effect of 
blurring the first order nature of the transition making 
measurements of Lq difficult. In order to reduce the ef¬ 
fects of tunneling we work with Lrms = \J{{L — Lq^). 
Due to fluctuations in L, Lrms will never reduce to zero 
even though (L — Lq) may become zero, however these 
fluctuations decrease as the lattice size increases. 

In order to measure the critical exponents, the relevant 
parameterization of the phase space must be identified. 


We therefore require an expression for the temperature 
T and external field H like directions in terms of /3p and 

h- 

In a magnetic system, a coexistence line occurs along 
the temperature axis. We therefore define the tempera¬ 
ture axis T for this system along the line of first order 
AC. 

Along this line the reduced temperature t is given by 


t = 




The definition of the external field H is less clear as 
there is no obvious preferred direction. In this work we 
used the definition used by Brezin and Drouffe in their 
mean field analysis of the critical exponents. 

The measurements of the critical exponents were car¬ 
ried out in two stages. Firstly we made a direct fit for L 
to the exponents. This method was applied for (3 and <5. 

Next we used finite size scaling 10 i which gave more 
reliable measurements of [3 and S and also provided a 
value for v and 7 . 


III. NUMERICAL RESULTS AND DISCUSSION 

1. Determination of P and v 


Figure ^ shows how Lrms changes with t for the differ¬ 
ent lattices. After discarding points with small t (where 
{Lrms) differs significantly from (L — Lq)), the best fit 
gives (3 = 0.41 ± 0.15. The solid line is the Ising value of 
/3 = 0.5. 



FIG. 2. Lrma versus the temperature 
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A more accurate measurement of (3 was obtained from 
finite size scaling. Figure ^ shows these results and com¬ 
pares the data with simulations run in the Ising limit, 
Pg ^ oo (solid line). The values for the critical expo¬ 
nents used in this graph are those of a four dimensional 
Ising model, /3 = 0.50 and ly = 0.50. Deviation from a 
single curve was observed ii v oi P were varied by more 
that 0.02. 



FIG. 3. Finite size scaling for p. L is the order parameter 
Lrms, I is the lattice size and t is the reduced temperatnre. 
The solid line represents the results from the Ising limit. 
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FIG. 4. Finite size scaling for 7 . y is the susceptibility, 
I and t are the lattice size and reduced temperature respec¬ 
tively. 


3. Determination of 5 

The measurement of 6 is less straight forward. It 
proved to be sensitive both to the location of the crit¬ 
ical point and to the lattice size (see Fig. |^). The main 
difhculty arises from the fact that L is not symmetric un¬ 
der the operation H —> —H generating different values 
either side of the temperature axis are = 4.4±0.4 for 
H > 0, and S = 3.0±0.3 for H < 0. Measurements of lat¬ 
tices with different sizes showed, that the shape of L vs.H 
becomes more symmetric as the lattice size increases. 


2. Determination of y 


Direct fits with 7 were unreliable and depended 
strongly on the dimensions of the lattices we used. 



Results from finite size scaling proved to be much more 
consistent as shown in Fig. Using the value of zz = 0.5 
found earlier gives 7 = I.O ± 0.1. 
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FIG. 5. Simulations passing through the critical point for 
lattice sizes 4“^ and 18^ 

The calculation of the range of 5 values took place sep¬ 
arately for each side around the zero of L. The combined 
result is (5 = 3.7±1.1. 



FIG. 6. Finite size scaling for 5 at (3p = 0.3194. L, I and 
t are dehned in figure 2. The curves are not symmetrical 
around the coexistence line. 

The asymmetry also affects finite size scaling. It 
proved impossible to collapse both H > 0 and H < 0 
data onto a single curve simultaneously as shown in Fig 
Treating each side independently gave v = 0.5 and 
(5+ = 2.5 and v = 0.5 and S~ = 4.7. 

4- Universality 

The measured values of the critical exponents show 
that the critical point in this model most likely behaves 
like a four-dimensional Ising model. The values of [3 and 
v fit the Ising values with very high accuracy, 7 has the 
Ising value too, but the measurement is less accurate. 

The only critical exponent that could be significant dif¬ 
ferent from the Ising value is 6. We ascribe our difficul¬ 
ties in measuring S to our definition of the external field 
direction H. It is not obvious from the model which di¬ 
rection is relevant and we were unable to isolate it from 
our data. The direction we adopted may therefore be 
a combination of T and H, thus our 5 measurement is 
probably contaminated with effects from (3. 


IV. SUMMARY 

This work is a numerically investigation of the behavior 
of the critical point in the Z 2 gauge lattice. The critical 
exponents were measured using Monte Carlo simulations, 
and were calculated directly as well as through finite size 
scaling. Other measurements were made in the Ising limit 
of the system. 

The calculated exponents had the four-dimensional 
Ising model values. The numerical measurements gave 
the same values as the theoretically calculated exponents 
of [3 and S. Two other critical exponents, v and 7, that 
were measured, gave also the Ising values. The finite 
size scaling for f3 created the same shape of curve as the 
one obtained from Ising data. The present high preci¬ 
sion runs set a high precision standard for the estimation 
of the character of other interesting critical points, no¬ 
tably the vortex/monopole condensation points in mixed 
action gauge theories and analogue points in supersym¬ 
metric theories. 
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